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ABSTRACT 

We consider the formation of OB associations from two perspectives: (a) the fractional gas con- 
sumption in star formation, e, per dynamical time scale tdyn in a galaxy, and (b) the origin of the 
so-called Kennicutt-Schmidt law that the rate of star formation per unit area is proportional to a 
power, a, of the surface density in H I and H2 gas when certain thresholds are crossed. The empirical 
findings that e 10^^ and a « 1.4 or 1.5 have simple explanations if the rate of star formation is mag- 
netically regulated. An empirical test of the ideas resides in an analysis of why giant OB associations 
are "strung out like pearls along the arms" of spiral galaxies. 

Subject headings: galaxies: ISM; galaxies: clusters: general; galaxies: magnetic fields; stars: formation 



I. INTRODUCTION 

Star formation in the low-redshift universe spans a 
huge range of linear sizes, from patterns on tens of kpc 
in spiral and irregular galaxies, to sub-pc cloud cores col- 
lapsing to a small fraction of an AU. Two empirical rules 
have brought integration to the subject: (a) the frac- 
tional consumption of gas per dynamical (e.g., rotation) 
time scale tdyn ~ 10® yr is a small number e ^ 10~^, and 
(b) the rate of star formation per unit area in galaxies is 

\ proportional to a nonlinear power of the surface density 
of cold gas in the galaxy: oc , with a « 1.4 (Martin 

. & Kennicutt 2001). 

One explanation of these two empirical relations 
invokes interstellar turbulence and stellar feedback 
(Elmegreen & Scalo 2004, Krumholz & McKee 2005). 
A problem is then the time required to damp turbulent 
fiuctuations, which is i/v^ where £ is the scale associ- 
ated with the velocity fluctuation v. This estimate is 
insensitive to the degree of magnetization of the medium 
(Stone, Ostriker, & Gammie 1998; Padoan & Nordlund 
1999). Prevailing theories postulate the turbulence to 
be virialized, which implies £/v ^ tg, where the free-fall 
time iff ~ tdyn ~ 10* yr for galactic distributions of cold 
interstellar gas on a scale of a kpc or so. If turbulence 
freely decays, the natural star-formation rate on galactic 
scales would be E* ~ Eg/tdyn, ~ 10^ times larger than 
the observed rate. Thus, the turbulence scenario faces 
two difficult obstacles, (a) pumping in energy at a suffi- 
cient local rate, usually by artificial means, to offset the 
natural turbulent dissipation (e.g., MacLow, Klessen, & 
Burkert 1998), and (b) explaining an efficiency of star 
formation e ^ 10~^ on local scales to offset the difficulty 
if 100% of the self-gravitating gas were to convert itself 
into stars on a time scale of tdyn ~ 10® yr. 

2. THE MAGNETIC APPROACH 

Fifty years ago Mestel & Spitzer (1956) noted that star 
formation would be impossible if the interstellar medium 
(ISM) were too heavily magnetized. The natural state of 
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the ISM would then give a rate of star formation which is 
zero. The debate between those attracted to magnetism 
(Mouschovias 1976; Nakano 1979; Shu, Adams, & Lizano 
1987; McKee 1989) and those favoring turbulence boils 
down then to whether it is easier to explain e ~ 1% 
starting from 0% or from 100%. 

Although the mass-to-fiux ratio does not have a well- 
defined meaning for a galactic disk of gaseous surface 
density E^ with an in-plane magnetic field of strength 
B, we introduce the concept of an ISM that is notionally 
subcritical in the competition between self-gravitational 
and magnetic forces: 



In the solar neighborhood. Eg w 5Mq pc^^, B k, b /iC, 
and we easily compute A„ ~ 0.3. Indeed, H I clouds in 
the Milky Way are all subcritical (Heiles 2004), whereas 
molecular clouds seem marginally critical (Crutcher & 
Troland 2006). 

Marginal sub-criticality of the ISM on average is prob- 
ably no coincidence (Shu et al. 1999). Heavily supercriti- 
cal gas clouds have long since collapsed to form stars. To 
induce the remaining clouds to collapse, any process has 
to overcome their notional stability. Successful processes 
only form stars with a low efficiency per unit collapse 
time scale because magnetic flux on a large scale cannot 
be destroyed but can only be redistributed. Thus, gas 
that becomes supercritical and collapses can only do so 
by making the remaining gas more subcritical and less 
prone to collapse. 

Under such conditions, the fraction that can collapse 
to form stars must depend on two time-scales: (a) tdyn of 
the large-scale dynamical processes that give rise to the 
proper circumstances for star formation, and (b) the age 
of the galaxy tgai. Their ratio gives the mean efficiency 
for star formation: e = tdyn/tgai ^ 10® yr/10^° yr ~ 1%. 
Galaxies with higher mean efficiencies in the past would 
have run out of star-forming clouds earlier. They have 
become galaxies with lower mean efficiencies today, not 
classified as spirals, bars, or magellanic-irregulars, but as 
ellipticals, lenticulars, or anemic spirals. 
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Unlike the notional A„, the real criterion for gravita- 
tional contraction on a large scale to form GMCs involves 
the surface density mn{2L) if we were to gather gas of 
mean volumetric number density n and mean molecu- 
lar mass m along a distance 2L in the direction of B 
(termed the "accumulation length" by Mestel 1985). In 
the absence of all forces other than self-gravity and elec- 
tromagnetism, this gathering produces a sheet, with the 
relevant component of -B|| then being the component par- 
allel to the sheet normal (Nakano & Nakamura 1978, 
Basu & Mouschovias 1994, Shu & Li 1997, Krasnopolsky 
& Gammie 2003): 



of motion as 



A; 



27rGi/2mn(2i) 



(2) 



T.g/2H, H is the 
B measured in 



With orientations as described, mn 
effective half-height of the gas, 
equation ([1]), and A = XnL/H. 

This line of thought (e.g., Lizano & Shu 1989) suggests 
that there are only two distinct ways by which gas that is 
notionally subcritical, A„ < 1, according to equation H]) 
can become actually supercritical, A > 1, according to 
equation Either decrease the denominator of equa- 
tion ([2]) by allowing the magnetic field gas to drift by 
ambipolar diffusion out of the lightly ionized gas in cold 
interstellar clouds (e.g., Nakano 1979), or increase the nu- 
merator of equation ([2|) by gathering the gas along mag- 
netic field lines over a parallel length 2L larger than the 
perpendicular dimension 2H to make A = X„L/H > 1. 

For the birth of massive stars in giant star clus- 
ters/associations occurring on dynamical time scales, we 
focus on the second process, in particular, on the mech- 
anism of transient gravitational instability in the post- 
shock regions of the arms of a spiral galaxy (e.g., Shetty 
& Ostriker (2006). We wish to extend the analysis to the 
implications for semi-empirical star-formation laws. 

3. COMPRESSION, CONTRACTION, COLLAPSE 

Consider a two-step process beginning with (a) a com- 
pression perpendicular to H and L across the field, such 
as by a magnetized spiral-density-wave shock (Math- 
ewson et al. 1972; Allen, Atherton, & Tilanus 1986). 
The compression (perpendicular to spiral arms) does not 
change the ratio of n compared to in a thin disk of 
height 2H, but it is followed by (b) a periodic contraction 
along the field (parallel to spiral arms), with expansion 
in between, which creates the effect of "feathering" in the 
gas distribution (see Kim & Ostriker 2006 and references 
therein). Dobbs & Bonnell (2006) also produce feath- 
ering in non-magnetic, non-self-gravitating calculations. 
The effect is ruined in pure hydrodynamic simulations 
when one includes the self-gravity of realistic amounts 
of interstellar gas (P. R. Woodward, personal commu- 
nication in the 1970s), which causes the collapse of en- 
tire spiral-arm segments (Chakrabarti, Laughlin, & Shu 
2003). Sufhcient magnetization can prevent this catas- 
trophe (Kim & Ostriker 2002). 

If the width of perpendicular compression (a) is W, 
the gathering of matter (b) along a length 2L in the 
parallel direction produces a mass M — 2LW'Eg. The 
magnetized flow through the arm makes the gravitational 
contraction along the length 2L a transient instability, a 
fact we temporarily ignore when we write the equation 



GM 



(3) 



The solution to equation ([3]) gives the time for complete 
contraction (into a feather) as 
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with no dependence on B because the contraction is 
along B. Given a post-compression example: = 
10 pc~2 (0.5 mag visual extinction for conventional 
dust-to-gas), W — 500 pc for the width of an arm, and 
2L = 800 pc for the half-spacing between giant H II com- 
plexes (at the head of feathers) strung along the outer 
spiral arms of some typical grand-design galaxies (M51, 
M83, MlOl), we get M = 4x lO^^M© and = 7x lO^r, 
a reasonable mass and age for a GMC complex (cf. Blitz 
et al. 2007). Contraction of L to a "final" value < H 
and W produces collapse in 3D. Surviving fluctuations 
( "turbulence" ) from the previous spiral-arm crossing in- 
troduces dumpiness that prevents actual complete col- 
lapse. 

li H — 100 pc, and A„ — 0.3, as in the solar neighbor- 
hood before and after perpendicular compression, then 
A = XnL/H = 1.2 after parallel contraction. Presum- 
ably, L is set by the time available locally to make A > 1 
marginally. Similarly, tc computed from observational 
data is comparable to the rotation time scale tdyn be- 
cause, apart from magnetism, the feathering instability 
has to overcome the tidal forces of the galactic gravita- 
tional field, the shear and rotation across spiral arms, 
and any turbulent velocities. This threshold is governed 
by some generalization of Toomre's (1964) Q parameter. 

With a resultant A = XnL/H greater than unity, the 
rate of GMC formation per unit area, EqmCj is given by 
^g/tc times the mass fraction ecMC of total gas converted 
into GMCs: 



^GMC = EGMC' 



tc 



4 [Gwy/\y, 



(5) 



Because half the gas is between feathers, we estimate 
eoMC 1/2- 

Being a mixture of supercritical cores and subcritical 
envelopes initiated by the turbulence of the recurring 
feathering instability, GMCs do not form stars with unit 
efficiency. The diverging mean flow in the interarm re- 
gion causes attached magnetic fields to tear the envelopes 
apart in the W direction faster than they can continue to 
accumulate mass in the L direction. To obtain the rate 
of star formation per unit area associated with GMC 
formation, we must multiply equation ([5]) by the mass 
fraction of gravitationally bound, supercritical cores that 
goes into continued gravitational collapse, ecoics, times 
the mass fraction of collapsing cores that results in the 
formation of stars, esp 1/3 (Alves, Lombardi, & Lada 
2007): 



E* — 6c 



sESpE 



GMC, 



(6) 

Empirically, Ccorcs ranges from 1% (Alves et al. 2007) 
to 10% or more (Alves, Lada, & Lada 1999); we adopt 
Scores ~ 6% as a representative value. Combining equa- 
tions ([5]) and ([6] gives a Kennicutt-Schmidt law, E, cx 
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Eg with a — 3/2 and a coefScient cx total efficiency 

^ ^GMC^corcs^SF 

estimated as (1/2) (6%) (1/3) = 1%. 
This estimate has large local scatter, but the argument 
of "natural selection" of §2 - that star formation on a 
galactic scale converts, in the mean, a significant frac- 
tion of the available gas into stars only on a time scale 
equal to the age of the system - ensures that e ~ 1% as 
a global average. 

For our conclusion on a to hold, W^^^ /L must have 
little dependence on T,g. The width W of gaseous spiral 
arms is set primarily by the properties of the background, 
stellar, spiral density-wave and not by (Roberts & 
Yuan 1970). The accumulation length 2L is set, we 
have already claimed, by the need to make the post- 
shock A marginally greater than 1. Observationally, 
L has much less variation from region to region than 
Sg. And the large-scale turbulence at ^ 8 km s~^, re- 
newed at each spiral-arm crossing, can have a decay time 



800 pc/8kms = 10** yr. 



The results ^ and ([6]) are independent of the sce- 
nario that motivated their derivation. Similar arguments 
should apply to the weak B-star formation that occurs 
in GALEX images of outer disks (Boissier et al. 2007) 
or to the"starbursts" in interacting galaxies like the An- 
tennae system, NGC 4038/9 (Zhang, Fall, & Whitmore 
2001). Nevertheless, the development of a rigorous the- 
ory of the process and observational tests of its predic- 
tions are most decisively carried out by a resurrection 
of the classical problem of why giant H II regions in a 
spiral or barred galaxy are "strung out like pearls along 
the arms" (Baade 1963). 



We are grateful to Art Wolfe for discussions and to 
Eve Ostriker for her illuminating papers with Kim and 
Shetty, which inspired this work. The questions raised 
by an expert referee helped greatly to sharpen the final 
arguments. 
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